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SECTION  I 


INTRODUCTION 


When  radio  waves  emanating  from  earth's  satellites  traverse  the  ionosphere,  they 
undergo  both  a rotation  of  the  plane  of  polarization,  i.e.,  Faraday  effect,  and  a Doppler 
frequency  shift.  Analytical  techniques  have  evolved  which  utilize  the  two  phenomena  for  the 
study  of  the  electron  content  in  the  ionosphere . 

Utilization  of  the  Faraday  method  which  necessitates  that  the  satellite-transmitted 
signal  be  linearly  polarized  has  been  demonstrated  by  Berlin  and  Papet-Lepine  (1970), 
Bertinetal.  (1966),  Blackband  (1960),  Checcacci  (1966),  Gol ton  and  Walker  (1971),  Kersley 
and  Taylor  (1974),  Klobuchar  et  al.  (1968),  Klobuchar  and  Whitney  (1966),  Lawrence  et  al. 
(1963),  Llszka  (1961,  1966),  Lyon  (1965,  1970),  Mendillo  et  al.  (1970),  Merrill  and 
Lawrence  (1969),  Munther  (1966),  Rao  (1967),  Roger  (1964),  Shmelovsky  et  al.  (1963),  Yeh 
and  Swenson  (1961),  and  Yuen  and  Roelofs  (1966). 

The  Doppler  method  which  requires  that  the  satellite  transmit  at  least  two  coherent 
harmonically-related  frequencies  has  been  employed  by  Bhonsle  (1966),  de  Mendon9a  (1962), 
Evans  and  Holt  (1973),  Millman  and  Anderson  (1968),  and  Ross  (1960a,  1960b). 

The  combination  of  the  Faraday  and  Doppler  methods,  often  referred  to  as  the  hybrid 
technique,  has  also  been  successfully  applied  to  ionospheric  electron  content  investigations 
by  Arendt  and  Soicher  (1969),  Burgess  (1963),  de  Mendon^a  and  Garriott  (1962),  and 
Golton  (1962). 

In  this  report,  an  evaluation  is  made  of  the  accuracy  of  several  of  the  analytical 
approaches  which  make  use  of  the  Faraday  and  Doppler  phenomena  for  determining  the 
ionospheric  electron  content.  This  is  accomplished  by  means  of  a simulator-computer  pro- 
gram which  consists  of  a satellite-orbit  generator,  a time-variant  three-dimensional  elec- 
tron density  model  and  an  earth  magnetic  field  model  expressed  in  terms  of  a series  of 
spherical  harmonics. 

In  Section  II  of  this  report,  the  analytical  formulations  for  deducing  the  electron  con- 
tent in  the  ionosphere  from  Faraday  and  Doppler  measurements  are  described. 
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A description  of  the  simulator  which  is  used  to  synthesize  the  Faraday  and  Doppler 
recordings  of  satellite  signals  is  given  in  Section  III. 

In  Section  IV,  the  analysis  of  the  simulated  Faraday  and  Doppler  data  and  the  accuracy 
results  of  the  various  methods  for  electron  content  determination  are  discussed. 

The  conclusions  of  this  study  are  presented  in  Section  V. 


SECTION  II 


THEORETICAL  CONSIDERATIONS 

2.1  FARADAY  ROTATION 

2.1.1  INTRODUCTION 

The  amount  of  angular  rotation,  n (in  radians),  experienced  by  a linearly  polarized 
wave  traversing  a one-way  path  in  the  ionosphere  can  be  represented  by  the  function 

h 

Ki  / 

Q = -Y  J H cos  6>  f(h)  Ng  dh  (2-1) 

^ 0 

4 

where  Kj  is  a constant  equal  to  2.362  x 10  cgs  units,  f Is  the  transmission  frequency  in 
Hz,  H is  the  magnetic  field  intensity  in  Gauss,  f(h)  is  the  secant  of  angle  between  the  ray 
path  and  the  zenith,  N^  is  the  electron  density  in  electrons/cm  , dh  is  the  height  differen- 
tial in  cm,  and  0 is  the  propagation  angle,  i.e.,  the  angle  between  the  direction  of  the 
earth's  magnetic  lines  of  force  and  the  direction  of  propagation. 

This  relationship  which  is  derived  in  Appendix  A contains  only  the  first-order  term 
for  the  refractive  index.  As  a first  approximation,  the  higher  order  terms  have  been 
neglected. 

Since  the  combined  function,  H cos  6 f(h),  varies  relatively  slowly  with  altitude 
(Millman  and  Rose,  1961),  it  is  valid  to  remove  the  terms  outside  the  integral  of  Equa- 
tion (2-1). 

Thus,  the  angular  rotation  can  be  written  in  the  form 
*^1  - 

~ M Nj  (2-2) 

where  llf  is  the  mean  value  of  the  geometric  magnetic  factor 
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and  is  the  integrated  electron  density  in  a vertical  column  up  to  the  satellite  altitude,  h^, 
h 

s 

f N dh  (1^—1) 

t J e 
0 

An  alternative  form  for  expressing  Equation  (2-3)  is 
K 

Q ^ -4  H cos  N (2-0) 

r 

where  N is  the  integrated  electron  density  along  the  ray  path,  i.e.,  oblique  path, 

R h h 

S 

= y dr  = J f(h)  dh  f(h)  j dh  (2-G) 

0 0 0 

For  any  ray  path,  the  parameter  M or  H cos  0 can  be  readily  specified.  In  this  study, 
the  magnetic  field  intensity,  H,  and  the  propagation,  angle,  0 , are  determined  by  assuming 
a spherical  harmonic  model  for  the  earth’s  magnetic  field.  A complete  description  of  the 
model  is  given  in  Appendix  A. 

According  to  Equations  (2-2)  and  (2-5),  the  electron  content  along  either  a vertical  or 
oblique  path  can  be  deduced  when  the  Faraday  polarization  rotation  angle,  is  known. 

In  general,  for  radio  waves  emitted  from  satellites  and  observed  on  the  ground,  the 
experimentally-measured  total  angular  rotation,  is  ambiguous  in  that 

= ± (n?r  i Af2)  (2-’?) 

where  n is  a positive  integer  and  Af?  is  the  acute  polarization  angle  which  would  normally 
be  indicated  in  a satellite-amplitude  measurement. 

2.1.2  SINGLE  FREQUENCY  METHOD 

It  is  possible,  however,  that,  for  transmissions  in  the  VHF  and  UHF  range,  the 
polarization  rotation  could  be  less  than  /2)  radians.  This  could  occur  under  certain  con- 
ditions; that  is,  when  observations  are  made  during  the  nighttime  or  at  certain  aziimilh- 
elevation  angle  orientations  toward  the  polar  ionosphere  (in  the  case  of  a polar  orl)iting 
satellite)  where  near  perpendicularity  with  the  earth's  magnetic  field  could  be  attained. 
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Assuming  that  the  angular  rotation  is  initially  less  than  (7i’/2)  radians,  then,  when  (tt/2) 
radian  rotation  does  occur,  the  integrated  electron  density,  i.e.,  electron  content,  can  be 
readily  derived  from  Equations  (2-2)  or  (2-5)  utilizing  a single  transmission  frequency. 

In  radar-lunar  studies  of  the  ionosphere  by  the  single  frequency  - Faraday  method, 
the  ambiguity  problem  was  resolved  by  theoretically  estimating  the  expected  angular  rota- 
tion along  various  earth-moon  paths  (Millman,  1964).  The  measured  angular  rotation,  as 
defined  by  Equation  (2->7),  that  best  correlated  with  the  theoretical  calculations  was  then 
selected  as  the  parameters  to  be  used  in  evaluating  the  integrated  electron  density  along  the 
ray  path.  In  the  theoretical  computation  of  the  magnitude  of  the  Faraday  rotation,  the  elec- 
tron density  distribution  derived  from  ionosonde  data  was  used  to  characterize  the  iono- 
sphere up  to  the  height  of  maximum  ionization  of  the  F-layer.  Above  the  peak  of  the  F-layer, 
the  distribution  of  electron  density  with  height  was  assumed  to  follow  a Chapman  model. 

2.1.3  TWO  FREQUENCY  METHOD 

One  method  for  resolving  the  ambiguity  of  the  complete  munber  of  polarization  rota- 
tions is  to  employ  two  closely-spaced  frequencies.  By  comparing  the  amplitude  fading  on, 
for  example,  40  and  41  MHz,  the  ambiguity  could  be  reduced  to  a multiple  of  20  half- 
rotations (Blackband,  1960).  Alternate  approaches  to  the  removal  of  the  n^r  ambiguity  in 
the  Faraday  rotation  measurements  have  been  proposed  by  Crooker  (1970),  Titheridge 
(1971),  and  Smith  (1971). 

When  two  harmonically- related  frequencies  are  used,  i.e.,  mfj^  = f^,  it  can  be  shown 
from  Equation  (2-2)  that,  when  neglecting  2nd  order  effects,  the  electron  content  Is  givenby 

(2-8) 

2.1.4  DIFFERENTIAL  POLARIZATION  ROTATION  ANGLE  METHOD 

The  ambiguity  problem  can  be  avoided  by  the  measurement  of  the  change  in  the  polar- 
ization angle  at  two  different  times,  t^  and  on  the  satellite  orbit  pass.  This  scheme 
which  is  often  referred  to  as  the  differential  Faraday  polarization  rotation  angle  method  re- 
quires that  the  integrated  electron  densities  along  the  two  different  ray  paths  be  identical. 
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This  assumption  implies  that  no  horizontal  gradients  of  electron  contert  exist  within  the 
time  interval.  Thus,  it  can  be  shown  from  Equation  (2-2)  that  the  electron  content  can  be 
obtained  from 
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2.1.5  POLARIZATION  ROTATION  RATE  METHOD 


(2-9) 


Another  method  which  can  be  used  to  determine  N^  is  one  which  involves  the  time  r..tv 
of  change  of  the  polarization  rotation  angle.  Differentiating  Equation  (2-2)  with  respect  to 
time  results  in 


n 


N^M 


(2-10) 


where  the  dot  signifies  the  time  derivative. 

Assuming  a horizontally  stratified  ionosphere,  then  N^.  = 0.  It  follows  that 
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A 
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(2-11) 


The  electron  content  can  therefore  be  determined  since  can  be  experimentally  measured 
and  M theoretically  predicted. 

A more  sophisticated  technique  which  eliminates  the  restriction  of  no  horizontal 
gradients  requires  the  2nd  derivative  of  which,  according  to  Equation  (2-9),  evaluates  to 
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If  it  is  assumed  that  N^^  is  of  the  form 
N^  = a t b t 


(2-12) 


(2-13) 


1 

I 

1 


I 


1 


2-4 


then  Nj  “ 0.  It  follows  from  Equations  (2-10)  and  (2-12)  that  the  electron  content  can  be 
expressed  by 


„ ^ ^ M n - 2M 
* ^1  MM-  2M^ 


2.1.6  LEAST  SQUARE  METHOD 

This  method  which  was  proposed  by  Burgess  (1963)  and  evaluated  by  Garriott  and 
de  Mendon9a  (1963)  involves  the  assumption  that  the  electron  content  can  be  represented  by 
a px)wer  series 

2 

Nj  = a+bt+ct  + ...  (2-15) 

At  t = 0,  which  corresponds  to  the  time  of  the  satellite  point  of  closest  approach, 

(N^)^  = a and  Equation  (2-2)  becomes 

<^t>o  = (2-16) 

where  A = (K^/f^). 

With  the  use  of  Equations  (2-15)  and  (2-18),  Equation  (2-2)  can  be  written  as 

n - 0^  = Aa  (M  - M^)  + Ab  Mt  + Ac  hft^  + . . . (2-17) 

It  is  noted  that  the  Faraday  rotation  angle  at  two  different  times,  n - 0 , can  be 

o 

readily  determined  from  the  experimental  data. 

The  least  square  error,  e , can  then  be  evaluated  from 


^ E[< 


(n-  - Aa(M-M^)-  AbMt  - AcMt  - 


whe’*e  n is  the  number  of  observations. 

Minimizing  < with  respect  to  the  coefficients,  a,  b,  c,  etc., 

9 ( _ _ - n 

9a  " 9b  " 9c  “ ” 


results  In  a set  of  linear  equations.  The  coefficients  can  then  be  evaluated  from  the  simul- 
taneous solution  of  the  linear  equations. 
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A method  of  solution  of  the  linear  equations  has  been  proposed  by  Garriott  and 
de  Mendonca  (1963). 

2.2  DIFFERENTIAL  PHASE  AND  DOPPLER  FREQUENCY  SHIFT 
2.2.1  INTRODUCTION 


The  differential  phase  or  dispersive  phase  method  can  be  considered  to  be  a modified 
form  of  the  Doppler  frequency  technique. 

The  differential  phase  between  two  harmonically-related  coherent  signals  transmitted 
from  a satellite  and  detected  on  the  ground  is  given  by 


_3 

where  is  a constant  equal  to  8.440  x 10  cgs  units,  and  a and  b are  constants  related  to 
the  satellite  transmitted  frequencies,  f ^ and  fg,  according  to 

afj  - bf^  (2-21) 


It  is  noted  that  differential  phase  between  the  two  frequencies  is  obtained  from 
A0  - a<t>^  ~ ^^2 


(2-22) 


where  0^  and  are  the  total  phase  shifts  encountered  by  the  radio  waves  traversing  a 
space-to-earth  propagation  path. 


Elquation  (2-20),  which  is  derived  in  Appendix  B,  is  based  on  the  assumption  that  the 
two  transmitted  signals  travel  along  the  same  ray  path  and  that  higher  order  terms  of  the 
refractive  index  can  be  neglected. 

The  differential  phase  can  also  be  written  as 


A0  = - 


f(h)  N^ 


(2-23) 


where  f(h)  is  the  mean  value  of  the  geometric  function,  f(h),  and  where,  as  indicated  by 
Equation  (2-6), 

Nj.  ~ f^  N^  (2-24) 
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The  absolute  phase  of  a satellite  signal  received  on  the  ground  cannot  be  measured. 
Thus  the  differential  phase,  Equation  (2-20),  the  relative  electron  content  along  the 

paths  from  the  satellite  to  a ground  terminal  throughout  the  satellite  pass  can  only  be 
inferred. 


2.2.2  DOPPLER  FREQUENCY  SHIFT  METHOD 

An  absolute  measure  of  the  electron  content  can  be  determined  from  the 
Doppler  frequency  shifts  of  the  coherent  transmissions  from  a satellite.  According  to 
Al'pert  (1958),  for  a satellite  moving  in  a circular  orbit  and  assuming  that  the  ionosphere 
is  nonrefractive  and  fixed  in  configuration,  both  in  time  and  space,  the  ionospheric  electron 
content  in  a vertical  column  can  be  evaluated  from 


(2-25) 


7 

where  Kg  is  a constant  equal  to  4.03  x 10  cgs  units,  h^  is  the  satellite  altitude  and  f^^^  and 
fj2  are  the  Doppler  frequency  shifts  of  the  satellite  emitted  frequencies,  f and  f^,  respec- 
tively. This  exp  ression  is  valid  for  satellites  located  at  high  elevation  angles  with  respect 
to  the  ground  receiving  station. 


As  shown  in  Appendix  B,  the  Doppler  frequency  shift  of  satellite  signals  traversing  the 
ionosphere  can  be  expressed  by 

K ^ 

/ "e'i'  <2-261 

0 

where  equals  1.343  x lO”^  cgs  units  and  R is  the  radial  velocity  of  the  satellite.  This 
equation  is  derived  on  the  basis  of  using  only  first-order  approximations  for  the  index  of 
refraction. 


Once  the  integrated  electron  density  is  known  for  one  point  on  the  satellite  trajectoiy, 
it  is  then  possible  to  convert  the  relative  electron  content  scale,  as  deduced  from  differen- 
tial phase  as  defined  by  Equation  (2-23),  into  an  absolute  measurement.  Thus,  the  electron 
content  in  a slant  column  can  be  derived  for  all  points  on  the  satellite  orbit  from 
Equation  (2-24). 
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2.2..!  DOPPI.KR  FRKQUKNCY  SLOPK  MKTIIOI) 

* idification  of  Al'pcrt's  expression  for  calculatintf  the  intenrateci  electron  density 
is  the  Doppler  slope  method  described  by  Arendt  et  al.  (19(55).  According  to  Arendt  et  al. 
(19G5),  when  the  satellite  is  at  the  point  of  closest  approach,  the  columnar  electron  content 
can  be  represented  by 


where  and  o-^  are  the  Doppler  slopes  at  the  inflection  points  for  the  frequencies  f ^ and 
respectively. 

It  should  be  noted  that  this  relationship  closely  resembles  Equation  (2-25)  derived  by 
Al'pert  (1958).  The  criterion  for  applying  the  Doppler  slope  method  has  been  discussed  by 
Arendt  (1966).  In  essence,  the  criterion  for  a minimum  prerequisite  for  the  applicability 
of  the  Doppler  method  is  given  by  the  inequality 


0 


k !2 

^d2  ^1 


< 1 


2.2.4  lONOSONDE  METHOD 


(2-28) 


The  ionosonde  method  which  was  suggested  by  Evans  and  Holt  (1973)  for  the  e.dibra- 
tion  of  the  differential  phase  records  utilizes  the  concept  that  the  electron  density  distribu- 
tion in  the  ionosphere  can  be  represented  by  the  Chapman  model  of  the  form 


N = N exp  - 
e m ^ 2 


(h  - h ) 

1 exp 


(h  - h ) 
' m 

H 


(2-29) 


where  H^  is  the  scale  height  of  the  neutral  particles  (atomic  oxygen)  and  is  the  electnm 

density  at  the  level  of  maximum  ionization,  h . 

m 


The  maximum  electron  density  of  the  layer  is  obtained  by  the  use  of  a vertical  inci- 
dence’ Jein^spherlc  sounder  and  is  related  to  the  ordinary  wave  - critical  frequency  of  the 
F-layer,  f^F2,  by 

""e  2 -8  ■> 

N = — ^ (f  f’2)  = 1.241x10  (f  F2)“  (2-30) 

m ^ ^ o o 

3 

where  N is  in  electrons/cm  , m is  the  electron  mass,  e is  the  electron  charge  and 
m e 

f F2  is  in  Hz. 
o 

The  scale  height  of  the  layer  (in  km)  is  defined  in  terms  of 

Hg  = 63.15  4 6.29  sin  +17.66  sin  ""  j (2-31) 

where  h is  the  local  time  in  hours  and  D is  the  day  of  the  year. 

This  analytical  model  was  derived  by  Klobuchar  and  Allen  (1970)  from  total  electron 
content  data  deduced  from  Faraday  rotation  measurements  of  radio  wave  transmissions  from 
the  ATS-3  geostationary  satellite.  The  observations  were  conducted  at  the  AFCRL  Saga- 
more Hill  Radio  Observatory,  Hamilton,  Massachusetts,  (geographic  coordinates:  42.6‘’N, 
70.8“W). 

In  order  to  employ  the  Chapman  distribution,  the  height  of  the  F-layer  maximum 
must  be  known.  From  incoherent  scatter  observations  made  at  the  Millstone  Hill  radar 
facility,  Westford,  Massachusetts,  (geographic  coordinates:  42. 6°N,  71.5°W),  Evans  and 
Holt  (1973)  have  concluded  that  h^  (in  km)  can  be  represented  by  the  function 

h^  = 280  4 40  cos  7T  j (2-32) 

When  the  three  parameters,  N , H and  h are  available,  it  is  then  possible  to  inte- 

m s m 

grate  Equation  (2-29)  to  obtain  the  integrated  electron  density  in  a vertical  column  to  the 
satellite  altitude.  Assuming  no  horizontal  gradients  of  electron  density,  the  electron  con- 
tent in  an  oblique  path  can  then  be  computed  from  Equation  (2-6). 


Wright  (1960)  has  shown  that,  if  the  electron  density  distribution  in  the  ionosphere 
is  of  the  Chapman  form,  the  total  Integrated  electron  density  can  be  expressed  in  terms  of 
the  scale  height  and  the  maximum  electron  density  by 


I N dh  = 4.133  N H (2-33) 

Jems 

0 

It  follows  from  Equation  (2-30)  that  the  electron  content  can  be  determined  from  the 
F-layer  critical  frequency  by 

N.  = 5.129  X 10®  (f  F2)^  H (2-34) 

t 'os 

It  should  be  noted  that,  according  to  Equations  (2-33)  and  (2-34),  the  maximum 
electron  density  or  the  critical  frequency  of  the  F-layer,  in  addition  to  the  scale  height, 
can  be  used  to  compute  the  electron  content.  However,  this  procedure  is  valid  for  satellites 
at  geostationary  altitudes  or  at  approximately  2 to  3 earth  radii  from  the  ground.  For 
TRANSIT  satellites  which  orbit  the  earth  in  the  vicinity  of  1000  km  altitude,  it  is  necessary 
to  integrate  Equation  (2-29)  to  the  satellite  altitude  for  electron  content  determination. 

2.2.5  LEAST  SQUARE  ME  THOD 

This  method,  which  is  similar  to  the  one  used  for  the  analysis  of  Faraday  rotation 
data,  was  also  proposed  by  Burgess  (1963). 

Expressing  the  electron  content,  N^,  as  a power  series  in  time.  Equation  (2-15),  then, 
at  t = 0,  which  corresponds  to  the  time  when  A<^  = 0,  (N^)q  = a and  Equation  (2-23)  becomes 


A0  = B f (h)  (NJ  = B f (h)  a 
o ''o'  to  ' 'o 

whe  re 


(2-35) 


(2-36) 
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Thus,  It  can  be  shown,  utilizing  Equations  (2-15),  (2-23),  and  (2-35),  that  the  differ- 


ential phase  at  two  different  times,  which  is  measurable,  can  be  written  as 


A0  - A<p^  = Ba  j^T(h)  - j ^ Bb  f(ir)t  + Be  T(H)  t^  + • • • 


(2-37) 


The  least  square  error  therefore  becomes 


n 

= ^ ^ j^(A0  - - Ba  ^)  - f^)  - BbT^h)  t - Be  f^)  t^  - • • • j (2-38) 


The  simultaneous  solution  of  the  set  of  linear  equations  obtained  where  < is  minimized, 
with  respect  to  a,  b,  c,  etc,,  will  result  in  an  evaluation  of  the  coefficients. 

A least  square  method  which  assumes  that  the  electron  content  varies  linearly  with 
time  has  been  considered  by  de  Mendonija  (1962). 

2.3  FARADAY  ROTATION  - DIFFERENTIAL  DOPPLER 


The  combination  of  the  Faraday  rotation  and  differential  Doppler  technique  to  deter- 
mine the  ionospheric  electron  content  was  originally  suggested  by  Burgess  (1962). 

The  mathematical  formulation  of  this  technique  can  be  readily  derived  by  differentia- 
ting Equations  (2-5)  and  (2-20)  with  respect  to  time.  Hence,  there  results 


n = — i r H cos  9 N + N H cos  6 I 

2 L r r J 


(2-39) 


■-?( 


■'a  fb‘  . 


Combining  the  two  expressions,  N^  can  be  readily  solved  from 


N = 
r 


H cos  6 


- 2tt  H cos  0 tr- 


72  2 

b - a 


(2-40) 


(2-41) 


where  Af^  is  the  differential  Doppler  frequency  shift 
Atd  - -i  Ai 


(2-42) 
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When  0,  Equation  (2-41)  simplifies  to 


(N  ‘ 

r o 


1 

H cos  0 


(2-43) 


It  should  be  noted  that  this  method  does  not  require  an  assumption  on  horizontal 
gradients  or  restricting  the  satellite  to  a circular  orbit,  i.e.,  no  vertical  velocity 
component. 


The  Burgess  hybrid 
de  Mendonca  and  Garriott 

I 


Faraday-Doppler  method  has  been  modified  by  Golton  (1962)  and 
(1962). 
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SECTION  III 

COMPUTER  SIMULATION  OF  SATELLITE  TRANSMISSIONS 

The  accuracy  of  the  various  analytical  techniques  for  determining  the  electron  content 
in  the  ionosphere  is  evaluated  by  means  of  a simulation-computer  program.  The  simula- 
tor basically  synthesizes  the  Faraday  and  Doppler  recordings  of  satellite  signals  that  would 
be  received  on  the  ground. 

The  logic  block  diagram  of  the  simulator  is  illustrated  in  Figure  3-1. 

In  the  satellite-orbit  generator,  the  azimuth  angle  (A),  elevation  angle  (E),  altitude 
(h),  range  (R)  and  range  rate  (R)  of  the  satellite  as  observed  as  a function  of  time  at  a 
ground-receiving  station  are  computed.  The  calculations  are  based  on  Kepler's  equations 
of  motion,  i.e.,  a two-body  orbit,  assuming  a rotating- spherical  earth  which  is  void  of  an 
atmosphere.  The  required  inputs  to  the  program  are  the  geographical  coordinates  and 
altitude  of  the  groimd  station  and  the  satellite-orbital  parameters  which  consist  of  the  three 
orientation  elements,  inclination  angle,  argument  of  perigee  and  longitude  of  ascending  node 
and  the  three  dimensional  elements,  semimajor  axis,  eccentricity  and  time  of  ascending  node . 

The  earth's  magnetic  field  is  represented  by  a series  of  spherical  harmonics.  In  this 
analysis,  the  set  of  80  spherical  harmonic  coefficients  for  Epoch  1965  derived  by  the  Inter- 
national Association  of  Geomagnetism  and  Aeronomy  (lAGA)  Commission  2 Working  Group  4, 
Analysis  of  the  Geomagnetic  Field  (1969),  is  used  to  specify  the  magnetic  potential  function. 
This  function  in  turn  1s  employed  in  the  computation  of  the  magnetic  field  Intensity,  H,  the 
inclination  angle,  I,  and  the  declination  angle,  D,  as  discussed  in  Appendix  A.  The  magnetic 
field  intensity  and  the  propagation  angle,  6 , which  is  a function  of  I and  D,  are  evaluated  at 
equal  increments  of  height  along  the  different  propagation  paths  to  the  orbiting  satellite. 

The  electron  density  profiles  generated  by  the  Penn  State  Mark  I Ionospheric  Model 
(Nisbet,  1970)  are  used  as  the  reference  ionosphere.  However,  the  simulator  CEin  also 
accept  as  an  input  any  spatial  distribution  of  electron  density.  The  required  li^)uts  to  the 
model  are  geographic  coordinates,  day  number,  time  of  day  and  solar  activity,  i.e.,  10.7 
cm  solar  flux  intensity.  The  output  of  the  model  is  in  the  form  of  hourly  values  of  the  elec- 
tron density  profile  ( Ng(h))  , between  120-  and  1260-km  altitude. 
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Figure  3-]  . I-ogic  IMock  Diagram  of  Faraday-Dopplcr  Simulator 


The  spatial  distribution  of  electron  density  is  attained  by  generatinK  vertical  profiles 
at  throe  locations  north  of  the  satellite-observing  site  and  three  locations  south  of  the  site, 
each  separated  by  10°  geographic  latitude.  The  data  for  the  six  additional  locations  yield 
latitudinal-electron  densitj'  coverage  and  information  on  the  existence  of  north-south  electron 
density  gradients.  Since  the  hourly  values  of  the  electron  density  profiles  can  be  converted 
directly  to  15°  longitudinal  values,  the  east-west  gradients  are  also  available. 


In  the  ray  path  - electron  density  generator,  the  electron  densities  along  the  different 
ray  paths,  [ N^(r)]  , are  determined  by  linear  interpolation  within  the  vertical  distribution  - 
electron  density  grid. 

In  the  Faraday  rotation  module,  the  first-order  term  of  the  polarization  rotation  angle 
(fi)  is  calculated  at  even  height  increments  along  the  ray  path.  The  output  data  also  consist 
of  the  geometric  magnetic  factor  (M),  the  magnetic  function  (H  cos  6 ),  the  integrated  elec- 
tron density  in  a vertical  column  (N.)  and  along  an  oblique  path  (N  ). 

V r 

The  differential  phase  between  the  two  coherent  transmitted  satellite  frequencies 
'v  (-^<h)  is  calculated  along  the  path  to  the  satellite.  In  addition,  the  Doppler  frequency 
shift  (f^|)  experience  by  the  two  frequencies  is  also  determined.  In  both  derivations,  the 
first-order  refractive  index  term  is  only  used. 

^ The  output  data  from  the  Faraday  rotation  and  differential  phase-Doppler  shift  modules 

which  in  reality  can  be  considered  to  be  a ground-satellite  receiving  system  is  processed  in 
the  analysis  program.  The  analytical  techniques  for  ionospheric  electron  content  determina- 
tion that  arc  evaluated  are  the  Faraday  single  frequency  method,  differential  polarization 
rotation  angle  method,  polarization  rotation  rate  method,  Doppler  frequency  shift  method, 
Doppler  frequency  slope  method,  ionosonde  method  and  the  hybrid  Faraday  rotation- 
differential  Doppler  method. 

The  electron  contents  deduced  by  the  analytical  techniques  are  compared  with  the  true- 
rcfcrcncc  values  obtained  by  integration  of  the  ray  paths  - electron  densities  in  the  Faraday 
rotation  module. 
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SECTION  I\' 

SIMULATOR  - DATA  ANALYSIS 

I 

'I'hc  simulated  Faraday  rotation,  differential  phase  and  Doppler  data  presented  in  this 
report  are  derived  on  the  premise  that  the  150-  and  400-MHz  transmissions  from  the 
: TRANSIT  satellite,  Object  No.  1970-067A,  are  recorded  at  the  General  Electric  Radio- 

^ Optical  Observatory,  located  at  42.  85°N  latitude  and  74,  07W  longitude  (54.  3°N  geomagnetic 

latitude)  near  Schenectady,  New  York.  The  satellite  which  was  launched  on  August  27,  1970 
is  in  a polar  orbit  with  an  inclination  of  approximately  90°,  apogee  of  1219  km,  perigee  of 
956  km  and  period  of  106.  9 minutes. 

The  earth  trace  of  the  satellite  orbit  at  approximately  1800  hours  GMT  on 
f June  3,  1974,  is  shown  in  Figure  4-1  together  with  the  geographic  coordinates  of  the  locations 

where  the  electron  density  - height  profiles  are  generated  utilizing  the  Penn  State  Iono- 
spheric model. 

The  altitude  and  azimuth-elevation  angles  of  the  satellite  orbit  as  predicted  for 
Schenectady,  New  York,  are  plotted  in  Figure  4-2.  The  azimuth  and  maximum  elevation 
angle  corresponding  to  the  point  of  closest  approach  is  276.  2°  and  71.  8°,  respectively.  The 
orbit  calculations  are  based  on  the  assumption  that  a vacuum  exists  between  the  ground  sta- 
tion and  the  satellite.  In  other  words,  tropospheric  and  ionospheric  refraction  and  time 
delay  effects  arc  neglected.  The  decrease  in  satellite  altitude  from  1200  to  1083  km  during 
the  oribital  pass  is  indicative  of  the  fact  that  the  satellite  has  a velocity  component  in  the 
vertical  direction.  In  the  evaluation  of  the  various  analytical  techniques  discussed  in  this 
report,  it  is  assumed  that  the  satellite  is  travelling  in  a circular  orbit,  i.  e. , vertical  velo- 
city component  is  neglected. 

The  vertical  distributions  of  electron  density  along  the  74.  07°W  meridian  at  1800  hours 
GMT  are  illustrated  in  Figure  4-3.  The  profiles  prertain  to  the  seven  locations  each  separ- 
ated by  10°  geographic  latitude.  The  center  profile  corresponds  to  the  receiver  site  at 
Schenectady,  New  York. 


i 

I 

i 

i 


4-1 


NORTH  LATITUDE  (DEGI 


i 

I 


It  is  evident  that  the  height  and  maximum  ionization  level  of  the  F-layer  peak  is  lati- 
tude dependent.  That  is,  the  electron  density  and  layer  height  are  a maximum  in  the  southern 
latitudes  and  a minimum  in  the  northern  latitudes.  In  the  vicinity  of  Schenectady,  the  F- 
laycr  peak  is  at  an  altitude  of  approximately  270  km. 

The  latitude  dependency  of  electron  density  signifies  the  presence  of  an  ionization 
gradient.  The  gradient  is  apparent  in  Figure  4-4  which  depicts  the  integrated  electron  den- 
sity in  a vertical  and  slant  (oblique)  column  to  the  satellite  altitude  and  along  the  satellite 
azimuth  - elevation  direction. 


The  simulated  Faraday  polarization  rotation  angles  of  the  150-  and  400-MHz  signals 
from  the  TRANSIT  satellite  are  plotted  in  Figure  4-5.  It  should  be  noted  that  the  satellite 
was  designed  to  transmit  circular  polarization  but,  on  a reception  at  the  ground,  the  polari- 
zation is  degraded  somewhat  to  elliptical  because  of  the  geometric  orientation  of  the  satellite 
antenna  with  respect  to  a ground  receiving  antenna.  For  this  analysis,  it  is  assumed  that  the 
TRANSIT  satellite  radiates  linear  polarization.  The  ambiguous  curves  are  representative 
of  the  data  that  would  be  experimentally  observed  while  the  unambiguous  curves  are  the  theo- 
retical estimates  of  the  total  rotation  that  the  two  frequencies  would  encounter  in  traversing 
the  ionosphere.  As  a result  of  magnetic  field  geometry,  minimum  angular  rotation  is  at- 
tained at  an  azimuth  angle  of  358.  2“  and  elevation  angle  of  9.  3°. 

In  evaluating  the  single-frequency  Faraday  rotation  method,  it  is  appropriate  to  con- 
sider only  the  400-MHz  data  in  order  to  insure  the  occurrence  of  minimum-angular  polari- 
zation rotation,  and  thus  avoid  the  ambiguity  problem. 


According  to  Figure  4-5,  the  ambiguous  angular  rotation  of  the  400-MHz  trans- 
mission is  less  than  tt  radians.  The  geometric  magnetic  factor,  M,  and  the  magentic 
function,  II  cos  0,  at  193.  5°  azimuth  and  9. 1°  elevation,  which  corresponds  to  the  direction 
of  the  (Jt/2)  radian  angular  rotation,  are  plotted  in  Figure  4-6. 


Utilizing  Equations  (2-2)  and  (2-5),  the  electron  content  along  a vertical  path,  .N^,  and 
along  a slant  path,  N^,  can  be  readily  deduced  for  a given  angular  rotation,  II. 


Figure  4-7  depicts  the  error  in  determining  the  vertical  and  the  slant  electron  content 
in  the  direction  of  ambiguous  Faraday  rotation  (7r/2)  radians  at  400  MHz.  The  error, 

AN^  j,  is  evaluated  from  the  expression 
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Figure  4-7, 


Error  in  Determining  the  Vertical  and  Slant  Electron  Content  in  the 
Direction  of  Ambiguous  Faraday  Rotation  of  tt/2  Radians  at  400  MHz 
By  the  Single  Frequency  Method 


1 3 

where  the  prime  signifies  the  theoretical  prediction.  For  this  example,  N'  4.  0221  x 10  ' 

,2  13  ^ 

electrons/cm“  and  N|  ■ 1.4050  x 10  electrons/cm“.  It  is  seen  that  zero  error  exists  for 

the  vertical  electron  content  when  the  value  of  the  M factor  at  :m  altitude  of  354  km  is  used. 

However,  for  the  slant  electron  content,  the  altitude  for  the  H cos  0 value  is  increased  to 

380  km.  It  is  noted  that,  in  the  (^/2)  radian  direction,  the  peak  of  the  maximum  ionization. 

is  located  at  approximately  280  km  altitude. 

In  the  case  of  the  vertical  electron  content,  the  error  is  increased  by  1 pei-ccnt  when 
the  M-value  altitude  is  varied  by  approximately  14  km.  For  the  shmt  electron  content,  when 
disregarding  the  reversal  portion  of  the  error  curve,  the  1 percent  error  change  is  associ- 
ated with  an  altitude  deviation  of  about  14.  5 km. 

In  determining  the  effect  of  the  magnetic  field  orientation  on  the  electron  content  - 
error  estimation,  the  electron  contents  were  computed  in  the  direction  of  the  (n  / })  radian  - 
ambiguous  angular  rotation  (20(1.5°  azimuth  and  44.  5°  elevation).  The  parameters.  M and 
H cos  n , in  this  direction  are  presented  in  Figure  4-8  as  a function  of  altitude.  It  is  appa- 
rent from  Figure  4-9  that,  for  both  the  vertical  and  slant  electron  content,  zero  error  is 
obtained  when  M and  II  cos  0 arc  evaluated  at  300  km  altitude  which  is  90  km  above  the 
maximum  ionization  level.  It  is  of  interest  to  note  that  the  altitude -error  slopes  of  the 

vertical  and  slant  electron  content  are  modified  to  10.6  km/1%  and  12.  1 km'T'7.  respec- 

13  ■’ 

lively.  The  calculations  in  Figure  4-9  are  based  on  1.  5086  x 10  electrons  'em"  and 

NJ  1.  1125  X 10^'^  clcctrons/cm^. 

For  general  analysis  of  Faraday  data  from  satellites  orbiting  at  altitudes  near  1000  km, 
Kersley  and  Taylor  (1974)  indicate  that  375  km  be  used  as  the  altitude  for  determining  M. 
However,  for  more  precise  imalysis,  they  recommend  that  the  height  should  !)e  so  km  aliove 
the  maximum  ionization  level. 

Titherifige  (1972)  is  of  the  opinion  that  a value  of  420  km  be  used  for  the  analysis  of 
geostationary  satcllite-Faraday  rotation  data  and  that  an  accuracy  of  <5  percent  in  electron 
content  estim.ation  can  be  expected.  Klobuchar  and  Allen  (1970)  on  the  other  hand,  assumed 
a constant  mean  ionospheric  height  of  350  km. 

•According  to  Yeh  (1974),  the  optimum  altitude  for  the  H cos  0 value  in  deducing  tlu- 
slant  electron  content  from  geostationary  satellites  is  680  km. 
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Figi  i-e  4-8.  The  Parameters  M and  H cos  6 in  the  Direction  of  Ambiguous 
Faraday  Rotation  of  tt/4  Radians  at  400  MHz 
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Jlant  Electron  Content  in  the 
>n  of  7t/4  Radians  at  100  MHz 


II. 


The  accuracy  in  measuring  the  polarization  rotation  angle  is  considered  to  be  on  the 
order  ±5°.  Thus,  for  an  angular  rotation  of  (Jr/4)  and  (Jr/2)  radians,  this  implies  an  addi- 
tional error  of  approximately  11. 1 and  5,6  percent,  respectively,  in  resolving  the  electron 
content. 

In  order  to  minimize  the  effects  of  ionization  gradients  on  the  accuracy  of  the  differen- 
tial polarization  rotation  angle  method,  it  is  preferable  to  employ  the  lower  frequency  to  de- 
duce the  electron  content  from  the  Faraday  data.  Since  the  Faraday  rotation  is  inversely 
proportional  to  the  frequency  squared,  as  indicated  by  Equation  (2-1),  a greater  number  of 
angular  rotations  occurs  at  the  lower  freqpency.  Thus,  for  a given  polarization  twist,  this 
permits  the  two  locations  on  the  satellite  orbital  pass,  as  specified  by  Equation  (2-9),  to  be 
spatially  closer  together. 

Although  Equation  (2-9)  expresses  the  relationship  for  determining  the  vertical  elec- 
tron content,  the  slant  electron  content  can  also  be  readily  derived  by  merely  replacing  the 
M factor  by  the  corresponding  H cos  6 terms. 

The  M and  H cos  0 parameters  along  the  ray  paths  in  the  direction  of  Faraday  rotation 

difference  of  {n  /2)  radians  at  150  MHz  are  plotted  in  Figure  4-10,  These  data  were  applied 

to  Equation  (2-9)  to  obtain  the  simulated  - computed  electron  contents.  Along  both  ray  paths. 

the  altitude  of  the  maximum  electron  density  was  located  at  approximately  270  km.  At  an 

azimuth  angle  of  346,  7°  and  elevation  angle  of  43.  0°,  the  true  vertical  (NJ)  and  slant  electron 

13  13  2 ^ 

content  (N')  were  1.0055  x 10  and  1.3939  x 10  electrons/cm  , respectively,  while  at 

^ 13  13 

227.  0°  azimuth  and  63.  0°  elevation,  they  were  1.  0679  x 10  and  1.  1825  x 10  electrons/ 

cm^,  respectively. 

In  deriving  the  errors  defined  by  Equation  (4-1)  and  shown  in  Figure  4-11,  the  mean 
value  of  N|  (NJ  1.  0367  x 10^^  electrons/cm^)  and  (N^  = 1.  2882  x 10^^  electrons/cm^) 
were  used.  The  interesting  disclosure  of  Figure  4-11  is  the  wide  separation  in  the  heights, 
252  km  and  613  km,  at  which  M and  H cos  9,  respectively,  are  evaluated  to  obtain  zero 
error.  The  slopes  of  the  error  curves  of  the  vertical  and  slant  electron  content  are  ap- 
proximately 6.  0 and  7. 1 km/1%,  respectively. 

In  order  to  facilitate  the  analysis  of  the  polarization  rotation  rate  method,  the  iono- 
sphere is  assumed  to  be  horizontally  stratified.  Equation  (2-11)  which  expresses  in 
terms  of  the  time  rate  of  change  of  the  polarization  angle  and  the  geometric  magnetic  factor 
can  also  be  employed  for  the  calculation  of  by  simply  substituting  ri  cos  9 for  M. 
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Figure  4-11.  Error  in  Determining  the  Vertical  and  Slant  Electron  Content 

at  150  MHz  by  the  Differential  Polarization  Rotation  Angle  Method 


The  geometric  magnetic  factor,  M.  and  the  magnetic  function,  H cos  i) . sliown  in 
Figures  4-12  ;md  4-13,  respectively,  illustrate  the  spatial  variation  of  the  parameters  at 
constant  altitudes  in  the  direction  of  the  satellite  trajectory. 

The  time  derivative  of  M and  H cos  0 plotted  in  Figure  4-14  and  the  magnitude  of  the 
terms  listed  in  Table  4-1  were  used,  in  conjunction  with  Equation  (2-11)  to  estimate  the 
accuracy  of  the  polarization  rotation  rate  method. 

As  shown  in  Figure  4-15,  the  altitude  for  computing  M and  H cos  n to  attain  zero 
percent  error  appears  to  be  a function  of  the  propagation  direction,  or  in  other  words,  the 
orientation  and  intensity  of  the  magnetic  field.  It  is  seen  that  the  mean  field  height  for  the 
slant  electron  content  is  on  the  order  of  300  to  400  km  above  that  for  the  vertical  electron 
content.  The  slopes  of  the  altitude  - vertical  electron  content  error  curves  arc  approxi- 
mately 1.  r>  and  fi.  3 km/1%  as  compared  to  the  corresponding  slant-slopes  of  4.  0 and  8.  0 
km/1%. 

The  simulated  relative  phase  difference  between  the  150-  and  400-MHz  transmissions 
received  on  the  ground  is  presented  in  Figure  4-16.  The  calculations  are  based  on  Equation 
(2-23)  utilizing  the  first  order  refractive  index  term  and  assuming  that  the  phase  comparison 
is  performed  at  150  MHz  with  the  constants  a 1 and  b (3/8). 

It  is  seen  that  the  differential  phase  undergoes  two  reversals.  The  first  reversal 
occurs  at  the  beginning  of  the  orbital  pass  in  the  high-altitude  ionosphere  while  the  second 
takes  place  at  an  azimuth  angle  of  287.5°  and  elevation  angle  of  71.  5°.  It  should  be  noted  that 
the  latter  appears  approximately  10  seconds  in  time  before  the  satellite  roaches  the  point  of 
closest  approach  which  is  located  at  276.  2°  azimuth  and  71.  8°  elevation.  The  time  difference 
between  the  two  data  points  is  attributed  to  the  presence  of  ionization  gradients. 

The  (Xjppler  frequency  shift  of  the  two  transmitted  signals  -^hown  in  Figure  4-17  are 
derived  on  the  basis  of  free  space  conditions.  Because  the  satellite  is  in  a noncireular 
orbit,  zero  Doppler  occurs  about  3 seconds  later  in  time  than  the  point  of  closest  approach. 
3'lie  zero  Doppler-angular  coordinates  arc  272.9°  azimuth  and  71.  77°  elevation. 

Figure  4-18  which  is  a plot  of  the  simulated  Doppler  frequency  shift  induced  hy  the 
ionosphere  was  computed  from  the  second  term  on  the  right  side  of  Equation  (2-26).  ,\s  in 

the  case  of  the  differential  phase,  the  higher  order  terms  in  the  refractive  index  ha\e  l)C('n 
neglected.  It  is  noted  that  the  spatial  position  of  the  zero-ionospheric  Doppler  frc<iu('ncy 
shift  coincides  with  the  second  phase  reversal  point  in  Figure  1-16. 
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Figure  4-14.  Time  Derivative  of  the  Parameters  M and 
H cos  6 
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TABLE  »-l 


PARAMETERS  USED  IN  THE  EVAl.l’ATION  OE  'niE  AUUUHAUY 
OF  THE  I>()LARIZATR)N  ROTATION  RATE  METHOD 


Parameters 


Azimuth  Angle  (Deg) 

Elevation  Angle  (Deg) 

150  MHz  Ambiguous  Polarization 
Rotation  Angle,  Q,  (Deg) 

d H/dt  (Rad/sec) 

rheoretical  V'ertical  Electron  Content 
(Electrons/cm^) 

Theoretical  Slant  Electron 
Content  (Electrons/cm2) 


A 

, H 

1. 1 

:155.  2 

314.  1 

22.  5 

07.  3 

43.  9 

41. 7 

0.  00«;i7 

0.  0071.5 

9,  7042  X 10^“ 

1.  0295  X 10 

2. 0244  X 10^^ 

1.  1002  X 10 

Figure  4>15.  Error  in  Determining  the  Vertical  and  Slant  Electron  Content 
at  ISO  MHz  by  the  Polarization  Rotation  Rate  Method 
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Figure  4-17,  Simulated  Free  Space  - Doppler  Frequency  Shift  of  150-  and  400-MHz  Transmissions 
From  TRANSIT  Satellite,  Object  No.  1970-067A,  June  3.  1974 


'riu'  conversion  of  the  relative  electron  content  obtained  from  differential  phase 
measurement  into  an  absolute  value  can  be  accomplished  by  means  of  the  Doppler  frequency 
shift  method  proposed  by  A 1 'pert  (1958), 

The  estimates  of  the  accuracy  of  this  method  for  four  satellite  locations  near  the 
point  of  closest  approach  are  contained  in  Table  4-2,  the  calculations  being  based  on 
Kquations  (2-25),  (2-26)  and  (4-1).  It  is  seen  that  the  error  in  determining  the  vertical 
electron  content  could  vary  from  -56.  5%  to  fll4.  4%.  It  is  possible  that  the  accuracy  de- 
gradation is  due  to  the  fact  that  the  satellite  is  not  in  a circular  orbit. 

'I'hc  Doppler  frequency  slope  method  of  Arendt  et  al.  (1965),  described  by  Equation 
(2-27),  can  also  be  applied  to  the  electron  content  calibration. 

Table  4-3  lists  the  values  of  the  parameters  at  the  point  of  closest  approach  used  in 
the  Doppler  slope  calculations.  It  is  found  that  there  is  a +29.  4%  error  in  the  estimation  of 
the  vertical  electron  content.  It  is  of  interest  to  note  that  100%  accuracy  would  be  attained 
if  the  satellite  is  assumed  to  be  at  an  altitude  of  889.  5 km. 

Vertical  incidence  ionospheric-sweep  frequency  soundings  in  conjunction  with  the 
Chapman  electron  density  profile.  Equation  (2-29),  can  be  employed  for  the  conversion  of 
the  differential  phase  data  into  an  absolute  electron  content  measurement  (Evans  and  Holt, 
1973).  In  appl3dng  the  Chapman  distribution,  it  is  necessary  to  infer  from  analytical  models 
the  scale  height,  H^,  and  the  height  of  the  F-layer  maximum  ionization,  h^.  Since  the 
simulated  receiver  site,  i.  e. , General  Electric  Radio-Optical  Observatory,  is  approximately 
211  km  west  of  the  Millstone  Hill  radar  facility,  it  is  valid  to  assume  the  models  for  H and 
h^.  Equations  (2-31)  and  (2-32),  suggested  by  Evans  and  Holt  (1973).  In  addition,  it  is 
assumed  that  the  ionosonde  can  accurately  measure  the  F-layer  ordinary  wave-critical  fre- 
quency, fyF2,  from  which,  according  to  Equation  (2-30),  the  maximum  electron  density  of 
the  F-layer,  N , can  be  deduced.  It  should  be  noted  that,  for  each  one  percent  error  in 

the  f F2  determination,  a 2%  error  will  be  introduced  in  N . 
o m 

The  results  of  the  analysis  of  the  ionosonde  method  are  summarized  in  Table  4-4. 

An  error  of  t^ll.  1%  implies,  that,  for  this  example,  the  ionosonde  method  would  give  a 
greater  electron  content  than  that  stipulated  by  the  Penn  State  model.  This  is  clearly  evident 
in  Figure  4-19  which  depicts  the  true  and  the  ionosonde -predicted  electron  density  profile. 


TAHLK  i-3 


KKUOU  KSTlMA'l'ION  OF  TIIK  OOl^Pl.FU  FltKCM  HNC Y 
SLOl’E  METHOD 


A/.imuth  Anglo  (Dog)  | 27(). « 


Elevation  Anglo  (Dog)  71.  h 

Sncollito  Altitudo  (Km)  1150.7 


150  MHz  Total  lX)pplor  Froquoncy  13.922 

Slopo  (Hz/soc) 

400  MHz  Total  lX)pplor  Froquonoy  37.132 

Slopo  (H/'Soo) 

1 3 

Estimatod  \’ortical  Electron  1.342H  x 10 

C ont out  ( E 1 00 1 rons /cm 2 ) 

Theoretical  \’ortical  Electron  1.03S0  x 10^'^ 

Content  (Electrons  cni“) 


Error  (Poroont) 
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KimOR  ESTIMATION  OF  THE  lONOSONDE  METHOD 


A/imuth  An^lc'  (l)t'g) 

27(i.  (i  : 

] 

! Elevation  Anglo  (Deg) 

71.8  ' 

j Satellite  Altitude  (Km) 

1 1 50.  7 ! 

1 

Local  Time  (Hours) 

1.4 

Day  Number 

154  j 

Seale  Height  (Km) 

8(1.  4 ' 

! Maximum  Ionization  Height  (Km) 

24.5.  4 

1 Maximum  Electron  Density 

4.  4 81  X 10''’ 

! (Electrons/cm'^) 

1 

.,14 

i Estimated  \ertical  Electron 

1.1,529  x 10 

('ontent  (Electrons  cm^) 

; 

1 

Theoretical  \ertical  Electron 

1. 0480  X 10 

Content  (Electrons/cm2) 

Error  (Percent) 

fl  1.  1 

I 

I 


I 
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LLtCFROIM  DENSITY  ( x 10  ELECTRONS/CM  I 


r’ifOJre  4-19.  Klectron  Density  Profile  at  Point  of  Closest  Approach 

of  TRANSIT  Satellite,  Object  No,  1970-007A,  1800  Iloui 
GMT,  -June  3,  1974  at  Schenectady,  New  York  as 
Predicted  by  the  Penn  State  Ionosphere  Model  and 
lonosonde  Method 


When  Faraday  rotation  and  differentia]  Doppler  signals  from  beacor  satellites  are 
simultaneously  recorded,  the  hybrid  method  proposed  by  Burgess  (1962)  can  be  applied  for 
ionospheric  electron  content  studies.  For  simplicity,  Equation  (2-43)  is  only  considered  in 
this  analysis.  This  relationship  is  applicable  for  the  position  on  the  satellite  orbit  at  which 
the  differential  IX)ppler  frequency  shift  or  the  time  rate  of  change  of  the  differential  phase 
is  zero.  It  should  be  noted  that,  on  replacing  H cos  0 with  M,  Equation  (2-43)  becomes 
identical  to  Equation  (2-11)  which  was  derived  for  the  polarization  rotation  rate  metliod 
assuming  a horizontally  stratified  ionosphere. 

Figure  4-20  is  a plot  of  the  time  derivative  of  the  parameters  M and  H cos  n in  the 
direction  of  the  differential  phase  reversal,  i.  e,  , 287.  6°  azimuth  and  71.  5°  azimuth.  The 
time  rate  of  change  of  the  polarization  rotation  angle  at  this  orientation  is  0.  00584  radian/ 
second. 

As  shown  in  Figure  4-21,  the  accuracy  of  the  hybrid  method  is  dependent  on  the  altitude 
at  which  the  mean  values  of  M and  H cos  0 are  computed.  Zero  error  is  achieved  for  the 
prediction  of  the  vertical  and  slant  electron  content  when  454  and  602  km,  respectively,  are 
selected  as  the  mean  field  heights.  The  vertical  and  slant  electron  content  error  slopes  are 
on  the  order  of  7.  0 and  7.  7 km/1%,  respectively,  which  are  similar  to  the  polarization  rota- 
tion angle  and  rate  results. 
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Figure  4-20.  Time  Derivative  of  the  Parameters  M and  H cos  0 
in  the  Direction  of  Differential  Phase  Reversal 
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Figure  4-21.  Error  In  Determining  the  Vertical  and  Slant  Electron  Content 
at  150  MHz  by  the  Faraday-Doppler  Hybrid  Method 


SECTION  V 


CONCLUSIONS 

It  is  demonstrated  that  the  concept  of  computer  simulation  can  be  applied  for  the  eval- 
uation of  the  various  anal3rtical  techniques  for  determining  the  electron  content  in  the  iono- 
sphere utilizing  beacon  satellite  transmissions. 

The  accuracies  of  the  Faraday  rotation  single  frequency,  angle  and  rate  methods  and 
the  Faraday-Doppler  hybrid  method  are  dependent  on  the  altitude  at  which  the  mean  geometric 
magnetic  factor,  M,  and  the  mear  magnetic  function,  H cos  0 are  evaluated. 

The  mean  field  heights  for  0%  error  in  the  estimation  of  the  vertical  and  slant  electron 
content  are  summarized  in  Table  5-1. 

The  deviation  in  the  M heights  with  respect  to  the  F-layer  maximum  ionization  level 
is  found  to  be  within  -50  and  +184  km. 

The  mean  field  heights  for  H cos  6 vary  from  approximately  150  to  425  km  above  that 
for  M except  in  the  case  of  the  Faraday  rotation  single  frequency  method  which  indicates  a 
height  variation  ranging  from  0 to  28  km. 

The  errors  in  predicting  the  electron  content  due  to  inaccuracies  in  the  mean  field 
height  are  listed  in  Table  5-2.  It  appears  that  a mean  field  height  error  has  the  least  effect 
on  the  accuracy  of  the  Faraday  rotation  single  frequency  method  and  the  greatest  on  the 
Faraday  rotation  rate  method. 

The  electron  content  prediction  accuracy  attained  by  employing  the  Doppler  phenomena 
is  summarized  in  Table  5-3.  Of  the  three  analytical  techniques  investigated,  the  ionosonde 
i method  seems  to  be  the  more  accurate. 

I 

I A more  thorough  investigation  is  required  to  determine  the  accuracy  that  can  be 

achieved  with  the  Faraday  and  Doppler  least  square  method  of  analysis  of  beacon  satellite 
signals. 
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TABLE  5-1 

ALTITUDE  OF  M AND  H COS  0 EVALUATION  FOR  ZERO  PERCENI' 
ERROR  IN  IONOSPHERIC  ELECTRON  CONTENT  ESTIMATION 


Analytical  Method 

Height  of  F-Layer 
Peak  (km) 

Mean  Field  Height  (km) 

M 

H cos  t) 

Faraday  rotation  single  frequency 

280 

354,300 

382.300 

Faraday  rotation  angle 

270 

252 

013 

Faraday  rotation  rate 

270 

220, 272 

540, 090 

Faraday-IX>ppler  hybrid 

270 

454 

002 

I 


i 

I 


TABLE  5-2 

ERROR  IN  IONOSPHERIC  ELECTRON  CONTENT  ESTIMATION  DUE  TO  MEAN 

FIELD  HEIGHT  INACCURACY 


Altitude-Error  Slope 
(km/1%  Error) 

Analytical  Method 

1 

Vertical  Electron 
Content 

Slant  Electron 
Content 

Faraday  rotation  single  frequency 

14.0,  10.6 

14.  5,  12.  1 

Faraday  rotation  angle 

6.  0 

7.  1 

Faraday  rotation  rate 

1.6,  6.3 

4.0,  8.0 

Faraday-Doppler  hybrid 

7.0 

7.  7 

TABLE  5-3 


i 


ERROR  IN  IONOSPHERIC  ELECTRON  CONTENT  ESTIMATION  UTILIZING 
DISPERSIVE  PHASE-DOPPLER  TECHNIQUE 


Analytical  Method 

Error 

(fo) 

Doppler  frequency  shift 

28,  8 - 114,  4 

Doppler  frequency  slope 

29,  4 

lonosonde 

11.  1 

r,-i 
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APPENDIX  A 


FARADAY  POLARIZATION  ROTATION  EFFECT 


When  a linearly-polarized  electromagnetic  wave  enters  the  ionosphere,  the  wave  sep- 
arates into  two  independent  components  both,  in  the  general  case,  elliptically  polarized  with 
opposite  senses  of  rotation.  For  frequencies  in  the  VHF  range  and  above,  the  two  compo- 
nents, i.e.,  the  ordinary  and  extraordinary  waves,  are  circularly  polarized. 

As  the  ionosphere  is  traversed,  the  two  waves  progress  with  different  velocities  of 
propagation  which  result  in  the  phase  relationship  between  them  to  be  continuously  changing. 
On  leaving  the  ionosphere,  the  circularly  polarized  components  recombine  to  form  a linearly 
polarized  wave  which  is  rotated  with  respect  to  the  original  linear  wave. 

The  differential  phase  shift  that  each  of  the  components  undergoes  in  an  element  of 
path  length,  dr,  is  given  by 

d0  = dr  (A-1) 

O A 

O 

and 

d0  = dr  (A-2) 

e A 

e 

where  the  subscripts,  oande,  refer  to  the  ordinary  and  extraordinary  wave,  respectively, 
and  and  A^  are  the  wavelengths  associated  with  the  phase  velocity,  and  of  the 
waves. 


By  definition 


and 


V 

po,  e 


o,  e 


(A-3) 


(A-4) 


where  c is  the  free  space  velocity,  f Is  the  transmission  frequency,  and  n^  ^ are  the  re- 
fractive indices  of  the  medium. 


Since  the  angular  rotation  of  the  plane  of  polarization  of  the  wave,  clf2,  as  it  travels 
the  disUmce,  dr.  is  given  by 

dQ  |(d0^-d0^)  (A-5) 

it  follows  from  Equations  (A-1)  through  (A-4)  that 

dfi  - ~ An  dr  (A-6) 


where  u.  - 27rf  and  An  = n - n . 

o e 

This  relationship  defines  the  angular  rotation  of  a linearly  polarized  for  a one-way 
propagation  path.  When  a two-way  transmission  path  is  considered,  an  additional  factor  of 
2 is  introduced. 


Thrj,  the  total  polarization  shift  (in  radians)  for  a one-way  propagation  path  in  terms 
of  a vertical  height  variable,  dh,  becomes 


Q 


Co 

2c 


An  f(h)  dh 


(A- 7) 


where  h and  h„  are  the  height  limits  of  the  path. 

X ^ 

The  function,  f(h),  which  is  the  secant  of  the  angle  between  the  ray  path  and  the  zenith, 
is  given  by 


f(h) 


r • h 
o 

^ i 

I h)^  - (r^  cos  E)^ 


(A-S) 


where  r is  the  radius  of  the  earth  and  E is  the  elevation  angle  of  the  antenna  beam, 
o 

According  to  Appendix  C,  the  difference  in  the  refractive  indices  between  the  ordinar>’ 
and  the  extraordinary  waves,  for  qua  si- longitudinal  propagation,  is 


An 


‘"N  ‘"H 


‘*'N  ‘*'H 


3 n 1 

COS  ^ ^ 2 


COS  0 


cos  0 
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(A-!)) 


where  is  the  angular  gyromagnetic  frequency  of  the  electron  about  the  earth's  magnetic 

field,  is  the  angular  plasma  frequency  of  the  ionosphere  and  0 is  the  propagation  angle. 

It  follows  that  the  one-way  polarization  shift  can  therefore  be  written  as 
^2 

n = — J'  cos  e f(h)  dh 

“n  ^ ^ f(h)  dh 

ojjj  cos  6 f(h)  dh  + • • • (A-10) 

2 

On  substituting  the  parameters  defining  and  from  Appendix  C,  this  relationship 
is  modified  to 


n = 


27rm  ^ c^  f^ 
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J"  Ng  H cos  0 f(h) 


dh 
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(A-11) 


where  e is  the  electron  charge,  m Is  the  electron  mass,  and  N Is  the  electron  density  and 

6 0 

H is  the  magnetic  field  Intensity. 


A-3 


When  utilizint;  the  numerical  values  of  the  constants  as  j;iven  in  Appendix  C,  th<' 
polari/.ation  rotation  in  cgs  units  is  simplified  to 


2.;i(>17  X 10‘^  r N H cos 
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cos  0 f(h)  dh 


1.H49:!  X 10 


17 


f 


O.SKiG  X 10 


11 


f 


H'^  cos^  0 f(h)  dh 


1 

h. 


/ 


cos  0 f(h)  dh 


(A-12) 


It  should  be  noted  that  the  first  term  in  Equation  (A-12)  is  generally  used  for  polari- 
zation rotation  calculations.  The  effect  of  neglecting  the  higher  order  terms  in  the  eom- 
plete  ex])ression  for  has  been  investigated  by  Yeh  (1960)  and  Ross  (1965). 

It  is  seen  in  Eciuations  (A-11)  and  (A-12)  that  the  magnitude  of  the'  polari/.ation  roUi- 
tion  angle  is  basically  a function  of  the  integrated  electron  density,  the  magnetic  field  in- 
tensity, M,  and  the'  ))ropagation  angle,  0 . 

It  can  be  shown  (Millman,  1959)  that  the  angle  0 is  given  by 

II  cos  ^ 1-cosf  sin  I - sin  f cos  I cos  (■;.  - D)]  (.\-i:i) 

where  I and  I)  are  the  magnetic  inclination  and  declination  angles,  respectively.  Tlu'se 
parameters  s|)ccify  the  direction  of  the  total  magnetic  intensity  vector  in  space. 

The  angle,  < , which  is  the  angle  between  the  ray  path  and  the  zenith  at  the  point  of 
magnetic  fit'ld  intersection,  is  defined  by 


A- I 


cos  K 


(A-U) 


It  is  noted  that  « is  related  to  the  function,  f(h).  Equation  (A-8),  by 


\L 


= sec  ^ f(h) 


(A- 15) 


The  angle  y which  is  the  geographic  azimuth  bearing  of  the  radar  location  as  measured 
at  the  subionospheric  point,  i.e.,  the  location  on  the  earth's  surface  directly  beneath  the 
magnetic  field  intersection  point,  can  be  derived  from  the  relationship  (Millman,  1958). 


y = tan 


sin  - ^p)  cos  <t>^ 


sin  cos  <pp  - cos  </)j^  sin  <pp  cos  - ^p) 


(A- 16) 


where  (p  and  \ are  the  geographic  latitude  and  east  longitude,  respectively.  The  subscripts, 
R and  P,  refer  to  the  radar  site  and  reflection  point,  respectively. 


The  parameters  I and  D can  be  deduced  from  a spherical  harmonic  model  of  the  earth's 
magnetic  field.  The  spherical  harmonic  representation  assumes  that  the  earth's  main  field, 
i.e.,  the  magnetic  field  that  excludes  such  phenomena  as  magnetic  disturbances  and  diurnal 
variations,  can  be  described  by  a regular  or  dipole  field  and  an  irregular  field. 


According  to  Chapman  and  Bartels  (1940),  the  magnetic  inclination  and  declination 
angles  for  the  spherical  harmonic  model  are  defined  by 


I = tan 


-1 


D = tan 
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(A-17) 
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where  X,  Y,  and  Z are  the  northward  horizontal,  the  eastward  horizontal,  and  the  downward 
vertical  component,  respectively,  of  the  total  magnetic  field  intensity,  H,  which  is  given  by 

i 

(A- 19) 
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The  components,  X,  Y and  Z,  are  related  to  the  magnetic  potential,  V,  by  the 
functions 
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(A- 29) 


where  r is  the  distance  from  the  center  of  the  earth,  0'  is  the  geographic  colatitude, 
p^'^'  (cos  0')  are  the  associated  Lengendre  functions  of  degree  n and  order  m,  and  g^^^^  and 
h^^^  are  the  coefficients  of  the  spherical  harmonic  expansion.  It  is  noted  that  the  first- 
degree  harmonic  terms  in  Equation  (A-23),  i.e.,  terms  with  n 1 and  m = 0,  reduce  to  th.al 
of  a dipole  potential. 
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APPENDIX  B 


DIFFERENTIAL  PHASE  AND  DOPPLER  FREQUENCY  SHIFT 


The  phase  of  an  RF  signal  received  on  the  ground  from  a transmitter  in  a satellite  can 
be  represented  by 


(^  = wt  - 


\ 


ti 
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where  <t>  is  the  phase  in  radians,  w is  the  transmitted  angular  frequency,  t is  the  time,  and 
X is  the  transmitted  wavelength.  The  parameter  P^.  is  the  phase  path  length  which  is 
given  by 
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where  n^^  and  n.  are  the  indices  of  refraction  of  the  troposphere  and  ionosphere,  respectively, 
and  dr  is  the  element  of  path  length. 

For  a satellite  which  transmits  two  harmonically  related  coherent  frequencies,  fj  and 
the  phase  of  the  received  signals  can  be  written  as 
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where  the  subscripts,  1 and  2,  refers  to  the  two  frequencies. 

The  nonionospheric  terms  can  be  eliminated  by  deriving  the  phase  difference  between 
the  two  signals  which  can  be  expressed  by 

A</)  = a 0j  - b 02  (B-5) 
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where  the  eonstants  a and  h are  related  to  the  frequency  ratio 
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Since  the  index  of  refraction  in  the  troposphere  is  independent  of  frequency  for  fre- 
quencies less  than  approximately  15  GHz,  the  phase  difference  evaluates  to 
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As  derived  in  Appendix  C,  the  ionospheric  refractive  index  is  given  by 
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where  is  the  plasma  frequency,  co^j  is  the  gyromagnetic  frequency  of  the  electrons 


II 


and  0 is  the  propagation  angle.  For  simplicity,  the  subscript  i is  no  longer  retained.  The 
positive  signs  of  the  ± terms  are  associated  with  the  ordinary  wave  while  the  negative  with 
the  extraordinai’y  wave. 
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Substituting  Equation  (B-8)  in  Equation  (B-7),  there  results 
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When  utilizing  Equation  (B-6),  this  expression  reduces  to 
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Further  simplification  can  be  attained  on  substituting  in  the  definitions  of  and 


ui„  given  in  Appendix  C.  It  follows  that 
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where  e is  the  electron  change,  m is  the  electron  mass,  N is  the  electron  density  and  H 

0 0 


is  the  magnetic  field  intensity. 

It  is  of  Interest  to  note  that  except  for  the  constants  a and  b,  the  2nd,  5th  and  last  term 
on  the  right  aide  of  Equation  (B-11)  are  identical  to  the  Faraday  rotation  terms  derived  in 
Equation  (A-11)  of  Appendix  A. 
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The  frequency  of  a radio  signal  emitted  from  a space  vehicle  and  received  on  the 
earth  experiences  an  apparent  shift.  The  phenomenon  often  being  referred  to  the  Doppler 
effect  occurs  because  of  the  relative  motion  between  the  transmission  source  and  the 
stationary  receiver  terminal. 

The  nonrelativistic  Doppler  frequency  shift,  f^,  in  the  ionosphere  can  be  defined  by 
the  relationship 
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On  substituting  in  the  definition  of  the  refractive  index,  Equation  (B-8),  there  results 
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It  can  be  readily  shown  that  Equation  (B-14)  simplifies  to 
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The  first  term  in  Equations  (B-14)  and  (B-15)  describes  the  Doppler  frequency  shift 
for  an  object  moving  in  free  space.  It  is  seen  that  the  Doppler  shift  is  a function  of  the  time 
derivative  of  the  integrated  electron  density  along  the  propagation  path  in  the  ionosphere,  the 
Faraday  polarization  rotation  and  the  refractive  bending  terms. 

The  magnitude  of  the  terms  in  the  Doppler  frequency  expression  has  been  evaluated  by 
Tucker  and  Fannin  (1968). 
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APPENDIX  C 
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IONOSPHERIC  INDEX  OF  REFRACTION 


The  index  of  refraction  in  the  ionosphere  in  the  absence  of  electron  collisions,  i . e . , 
absorption,  is  defined  by  Ratcliffe  (1959) 
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where  0 is  the  propagation  angle,  i.e.,  the  angle  between  the  magnetic  field  vector  and  the 
direction  of  propagation. 

The  parameter,  is  the  angular  gyromagnetlc  frequency  of  the  electron  about  the 
earth's  magnetic  field  and  is  given  by 
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where  H is  the  magnetic  field  intensity  in  gauss,  e is  the  electron  charge  (4,8  x 10*^®  esu), 

—28  10 
IS  the  electron  mass  (9. 1 x 10  gm)  and  c is  the  free  space  velocity  (3  x 10  cm/s). 

The  term,  is  the  angular  plasma  frequency  of  the  ionosphere  and  Is  given  by 
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where  is  the  electron  density  in  electrons/cm  . 
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It  should  be  noted  that  there  are  two  values  for  the  refractive  index,  Equation  (C-1). 
The  jxisitive  sign  is  associated  with  the  ordinary  wave  while  the  negative  sign  with  the 
extraordinary’  wave. 

The  quasi-longitudinal  mode  of  propagation  can  be  represented  by  the  condition 


It  can  be  readily  shown  that,  in  order  to  maintain  about  a factor  of  50  between  the  two 
terms,  the  approximation  is  valid  when  the  propagation  angle  is  restrained  between 
0°  _<  86.0°  at  100  MHz  and  0°_<  6»  _<  89.6°  at  400  MHz. 

For  the  quasi-longitudinal  case,  the  refractive  index  simplifies  to 
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Substituting  Equations  (C-2)  and  (C-4)  in  Equation  (C-8),  it  follows  that 
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The  condition  for  quasi-transverse  propagation  is  denoted  by  the  inequality 
(1  - X)^ 
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which  is  merely  the  reverse  of  Equation  (C-7).  For  this  case.  Equation  (C-1)  reduces  to 
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where  the  subscripts,  o and  e,  signify  the  ordinary  and  extraordinary  waves,  respectively. 


Referring  to  Equation  (C-8),  the  refractive  index  for  the  quasi-longitudinal  mode  of 
propagation  can  be  written  as 

n^  = 1 - X(1  ± (C-1 

Expanding  by  the  binomial  theorem,  there  results 


Further  expansion  can  be  accomplished  by  letting 
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In  order  words,  the  index  of  the  refraction  becomes 
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This  expression  can  also  be  written  in  terms  of  the  physical  constants  for  X and  V 
Kquations  (C-2)  and  (C-4).  It  follows  that 
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The  difference  between  the  refractive  indices  of  the  ordinary  and  extraordinaiy  waves 
is  therefore 
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